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Abstract 


Single series forecasting (SSF) is the forecasting of a 
Single time series based on the current and past values of 
the series. only. It 1s implemented for process control 
using a process model, a forecaster and a controller. The 
difference, or residual r(k), between the actual process 
output y(k) and the model output ¥(k) is fed back through 
the forecaster to the controller. Based on the present and 
past values of r(k) only, the forecaster generates forecasts 
r(k+A), where 2d is not less than the process plus sampling 
delays. Any controller may be used. However, Since the 
process model is already available, a predictive controller 
is recommended because it generates the appropriate control 
action at time k to compensate for the forecasted error at 


time k+). 


SSF was investigated as a means: 

1) to improve the performance of both conventional and 
modern control techniques, and 

2) to provide insight into the structure and 
performance of techniques such as Internal Model 


Control (IMC) as proposed by Garcia and Morari. 


It has been shown that, given perfect forecasting, SSF 
can result in perfect control even with the presence of 
unmeasured disturbances. If the available process model is 
Not Ppertect}e thenrpsshecheips Glomacompensate for both 


unmeasured disturbances and modelling errors. 
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The quality of control using SSF depends on the quality 
of the forecast and hence deteriorates in the presence of 
large amounts of random noise and/or strong nonlinearities. 
Furthermore, a perfect predictive controller will compensate 
for the effect of a disturbance in A time intervals when 
equals the sampling plus process delays. Thus, the best 
that SSF can do when combined with a perfect predictive 
controller is to achieve compensation at time k rather’ than 
kgs. The practical JUStVoacataon for actual SSF 


applications is therefore limited. 


Comparison of the SSF approach with other classical and 
modern control techniques does, however, provide some 
important insights. 

1) Classical proportional plus derivative feedback control 
can be interpreted as proportional control based on SSF 
of the control error or process output. 

2) IMC can be interpreted as an SSF system that uses the 
current value of the residual r(k) as an estimate of the 
future value r(k+dA). Futhermore, the design rule that 
the IMC controller should approach the inverse of the 
process model is readily seen as being equivalent to the 
classical design observation that perfect control is 
approached as the feedback gain approaches infinity. 
(There are of course stability and performance 
limitations with most practical systems.) 

3) IMC can be interpreted as a classical feedforward scheme 


based on the forecasted (estimated) disturbance rather 
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than the actual measured value. 


Simulations results also showed that SSF could improve 
the performance of adaptive control techniques, such as’ the 
Adaptive Predictive Control System (APCS), which use the 
magnitude of the estimation error as a criterion for 
Switching identification on and off and have no explicit 
identification of the noise term. During the periods that 
the adaptive parameter estimation is turned off, SSF is used 
to model the structure of the residual (disturbances plus 
modelling error) and augment the predictive control action. 
This second level of control can be justified on the basis 
that it prevents unnecessary corruption of the process I/0 
model by unmeasured disturbances and can reduce the control 


error. 
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Nomenclature 


Alphabetic 


a ARMA model parameters 

b ARMA model parameters 

C Feedback controller 

d Disturbances 

a Residual or estimated disturbances 
F Filter 

Gc Controller 

Gd Disturbance transfer function 
Gm Process model 

Gp Process transfer function 

Gx Feedforward controller 

k Sampling instant 

Ke Proportional constant 

r Residual 


Estimated residual 


FY 


R Reference model 

u Process input 

x Process input/output vector 

y Outpue 

y Estimated output 

Ya Desired output 

Vin Estimated output 

vs Setpoint 

va Disturbance effect on the output 
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Greek 


a Coefficient of AR forecaster 

Aa Upper bound on the absolute value of A(k) 
A Minimum upper bound of A, 

€ Control error 

6 Process parameters 

0 Model parameters 

r Process time delay plus sampling delay 
é Disturbances 

Ta Derivative time constant 

Ti Integral time constant 

} Process input/output vector 

w Switching factor 


Xvi 


2tesemaxrsa eee ne019 
/ > % 
, lassijeaigisq Lebom 


welsh pniiagice aula yesh emis aesz079 


Jrazenes eMiza svisevised 


trpiates sat? Leigeint 


{decoy suctue Mon? aneto7s 


= 


Toraes ontaaared 


sa 


1. Introduction 


1.1 Introduction 

Inkehempeal@ process Scontrol),S@process? lagsslandletime 
delays present’. a challenging: problem *to SmostY*’ control 
techniques. Conventional feedback control is still the most 
popular control technique despite the vast development in 
modern control techniques, but it performs poorly for 
systems with large time delays. 

Predictive control techniques have been developed to 
address this probelm. An example of an early development is 
the Smith Predictor. Modern control techniques such as’ the 
Self Tuning Regulator (STR), Adaptive Predictive Control 
System (APCS) and Internal Model Control (IMC) are all 
predictive in nature. Control action is determined based on 
the prediction of the process output by an input/output 
(I/O) model. 

It is obvious that the use of a perfect I/O model for 
prediction can improve control in the presence of process 
lags and time delays. However, non-idealities such as 
modelling errors and unmeasured disturbances generally 
degrade the performance. Unmeasured disturbances cannot be 
detected until their effects appear in the output of the 
process. In order to compensate for unmeasured disturbances 
before they affect theS processus output;e somesriorm: of 
forecasting of their effects on the output is necessary. An 


estimate of the current disturbance effect can be obtained 


ary 
oe an J 
sain tat 
omit bis epel. s8es07g eee: éspcoxg Isoimeda al 
oegeatq «6 eyaleb ” 


fovanes deem.) oF 
$aom set firte ef Los aod aageioal Lend Inevaes /eoupindoed 

Pees Poe ; 
nt snahyoleveb tesv sd? s3iq@eeS supiatice: -foxtnon asiuqog 


rs 2 ee e ; —_ 
i. jud. \eelpindses fornaos nsebou 


102 -yizocgd. .amvomtsg 7 
_ Joe's 6b sait epie: Agiw exsseya 
+ a _ 
63 Begofeaveh reed soved e5ucindosd lessi6e evistaties9 


“i 


o} 4heamol a: cb vfs as Qo siamaye RA’ .miederq eit eaerbba 7 


lovsas2-evi sorbet av Preses ated naidlaicuein enimuT ile - 
Ife 4%6 (ol) Logaied ieboM lenvegnt Snes) teaHA mateye 
1G bse6d 59 titel sO" ee, Gagine idashean a toten fi evizoibeq: 
cuasuC \suqndi Es ve 2003u0 saes014 Si3 To foiijntber39 odd 
. febom i 
a7 Lepon Oa) #293398q wide e2y 4hd Gand avoivde at 3 
efeue2g io acnssesm edd nic fosico>” Sve=qmc fee nobtoibexg 


7 


svevoH  ja¥eivh amit bas wpa! 


Hh 
A 
, 
) 
“ 
»- 
te 
ap 
Sy 
bas 
I 
a 
~ 
=! 
af 


ylikssney escead nang ee bs riegaiant bie: etotis - peas 
ga tonnss 286 cusdsuae iy! baieReanl ’. (Sngmye aa a 
sits In 3 jugs mo. SH a ni nesags. esate: sits tian & 
moneda tb ti ence 104. sdsensqines oF 198 ee 


_ 


ges — ne 
do od 699% 
cle 


from the difference between the measured process output and 
Che @eprediicted@soutput priromy theeer/0 ismodel /itBased ontthe 
current and past values of this estimate, future values of 
the disturbance can be obtained using single series 
Forecasting (SSH )RIE.ThesworkSingrhismethesise useudevotedmato 
investigating this use of SSF for compensating unmeasured 


disturbances. 


1.2 Objectives of the Thesis 

The overall objective of the thesis is to investigate 
the use of SSF in process control. The original interest of 
this project was to improve the control performance of 
adaptive control techniques by SSF. The scope was later 
expanded to include Internal Model Control (IMC) and 
conventional feedback and feedforward control. Conventional 
control techniques have been included because some of their 
fundamental principles can be used to interpret IMC and to 
introduce the use of SSF for process control. 

The specific objectives of this thesis are: 

1) to investigate the use of. SSF aS a means to improve the 
performance of both conventional and modern control 
techniques, and 

Pie toe gainSfinsighteeintotthe structure andtper formance of 
techniques such as Internal Modern Control (IMC) as 


proposed by Garcia and Morari (1982). 
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Since time series forecasting is a fairly well 
established area, the primary concern in the first objective 
is the viability of using SSF with existing control 
techniques rather than the relative merits of different 
forecasting methods. Perfect forecasting is assumed in many 
cases to establish the best improvement obtainable from SSF. 
Then a linear forecaster or an autoregressive forecaster is 


used to establish a more realistic case. 


1.3 Structure of the Thesis 

The contents of the seven chapters of this thesis have 
been arranged to give the readers a logical view of the 
development of the work on SSF. Following the introduction 
of the objectives of the thesis in Chapter One, Chapter Two 
explains the concept of SSF and presents some common 
forecasting methods. Readers who have prior exposure to SSF 
will be able to omit Chapter Two or use it aS a review. 

The main results of the thesis are presented in 
Chapters Three to Five according to the various control 
techniques which have been studied. Fors each control 
technique, a description of the technique is given followed 
by an analysis or an interpretation of the technique. Then 
the use of SSF with the control technique or the relation of 
SSF with it is examined. Finally, some numerical examples 


are given to illustrate the ideas presented. 
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The control structure for implementing SSF is 
introduced from a feedforward control point of view in 
Chapter Three. This structure serves as a reference for the 
use of SSF with IMC and APCS in Chapters Four and Five 
respectively. Also in Chapter Three, some fundamental 
principles of conventional feedback and feedforward control 
are presented explicitly as concepts for the purpose of easy 
reference by Chapter Four in the interpretation of the IMC 
structure. The thesis is deliberately structured to 
emphasize that many of the features of "advanced" control 
schemes such as IMC follow directly from a sound 
understanding of classical feedback and feedforward control. 

Chapter Six contains an application of SSF on a 
Simulated double effect evporator as a realistic evaluation 
of SSF. Chapter Seven presents conclusions from the results 
in the previous chapters and recommends ideas for future 


work. 
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2. Single Series Forecasting (SSF) 


2.1 Introduction 

The technique of single series forecasting (SSF) has 
been used extensively in economics, business, engineering 
and many other areas. The purpose of this chapter is to 
explain briefly the concept of SSF and present some simple 
forecasting methods. Special emphasis has been placed on 
the differences between SSF and prediction using an 
input/output (1/0) model. This chapter is intended for 


readers with no prior knowledge of SSF. 


2.2 Concept of SSF 

SSF is the forecasting of a single time series based on 
knowledge of the current and past values of the series only. 
To illustrate the idea of SSF, consider a meaSured discrete 
time variable r(k) as plotted in Figure 2.1. At a certain 
Sampling instant k, it is required to forecast the future 
values f of r. Assuming the factors affecting r are either 
unknown or unmeasured, the only information available for 
the forecast are the current and past values of Ce 
Therefore, estimates of the future values, f, have to be 


generated by SSF. 
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Variable r(k) 


Measured | Forecasted 
Values | Values 
k Time 
Figure 2.1 Illustration of Single Series Forecasting 


The term "forecasting" is used here to distinguish SSF 
from POnedicevon s uSingesan 1/0 model. An 1/0 model 
represents a functional relationship between two variables - 
EhCemInDUL son essthe = OutpDUL. If the model is accurate and 
there is no other factor affecting the output, perfect 
prediction can be obtained knowing the current and past 
input and output values. Therefore, SSF and I/O model 
prediction are different in terms of the information they 
USentoOLresprediction. 

To use SSF, it is obvious that the time series must 
have some deterministic structure Since it is impossible to 


accurately forecast a random variable. The quality of the 
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forecast depends on the dynamics of the time series, the 
lead time of the forecast, and the type of forecasting 
method used. The type of forecasting method used will be 
determined by the requirements of the individual 
application. Some common forecasting methods will be 


presented in the following section. 


2.3 Forecasting Methods 

The problem to be solved here is to forecast future 
values of a time series using the current and past values of 
the series only. Many methods based on different approaches 
are available. Three commonly used methods will be 


presented in the this section. 


1) Forecasting using current value 
There is actually no real forecasting done in this 
method because it uses the current measured value as the 


forecasted value. For a one-step ahead forecast, 


ECKtdy =er(k) a2ia9 


This method is satisfactory if the time series changes 


Slowly. It can also be used aS a conservative approach for 


time series containing a large amount of noise. 
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2) Linear extrapolation 
This method makes use of the first derivative of the 
time series at the current value to calculate the future 


values. The equation for a one-step ahead forecast is 


dr 
r(k+1) = r(k) + — At Wee 
dt 
where At = sampling time 


The derivative can be calculated by the Newton's 
backward difference formula. For example, the two-point 


formula and the three-point formula are given as follows. 


dr cék) a=En( k= 1) 


Two points; — = ——————— V2e3a) 
dt At 
dr 3r(k) - 4r(k-1) + r(k-2) 
Three points: | — = ————___—_ (2743p) 
at 2At 
In general, the use of two or three points is 


sufficient for calculating the first derivative. 

The advantage of this method is that the computation is 
Simple and the extrapolation is accurate for time series 
without sharp changes in the values. The disadvantage is 
that it is easily confused by noise. It is also not good 
FOL wlong termeLorecastingd. sin) general; sonly  sLonecasts ror 


one or two sampling periods ahead are reliable. 
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3) Forecasting using an auto-regressive equation 
In this method, a general aA-step ahead forecast is 


given by the auto-regressive (AR) equation. 


P(KtA) =taqrtk) + aer(k-1) + ... + a,r(k-ntl) C24) 


Phem@ coetficientS@,.a;, a2, ..3., @, May or may not be 
known. If they are unknown, they have to be identified 
either off-line or on-line. Then the equation can be used 
to calculate the forecasted values. 

To use the AR equation for disturbances whose nature 
generally changes with time, the coefficients have to _ be 
identified on-line by a parameter identification algorithm. 
If the time series has a structure, the AR equation can be 
interpreted as a model for the series. The coefficients 
will get closer to the true values as the identification 
process continues. If the coefficients converge to the 
correct values, the forecast will be perfect. However, it 
is possible to obtain accurate forecast without using 
accurate coefficients. The primary objective is to have 
accurate forecasting rather than accurate coefficients. 

There are many on-line parameter identification 
algorithms available. Two of them will be described here 


because of their simplicity and fast convergence. 
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a. Recursive least square algorithm 
The time series of the variable r is assumed to be 
described by the equation 


r(k) }(k-1) @ (225) 


where 
Gee= (Sa7 (6 Goo one | 


SCk=1) eae leer (kod) r(k-2) ee pe riken) 3] 


The recursive least square algorithm, e.g. Astrom and 


Eykhoff (1971), is given by 


Oe oO Cko ete Kh Ki rT Gky® Ck 1).60k-3)3) (25) 


K(k-1) P(k-1)@*t(k-1) [1+®(k-1)P(k-1)8*(k-1)]-' (2 7,) 


1 
PC kin ft P Oka) —K (k=) 
[1+@(k-1)P(k-1)*t(k-1) ]K*(k-1)} (258) 


6B is the exponential forgetting factor whose value is 


chosen as 0.9 < B < 1. P(k) is the covariance matrix. 


b. APCS projection algorithm 

The Adaptive Predictive Control System (APCS) was 
proposed by Martin-Sanchez (Aa. The recursive 
identification algorithm, also known as a projection or a 


gradient search algorithm, used by APCS is 
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identification on and off. A more detailed definition of 
W(k) will be given in Chapter Five. 

Forecasting with an AR equation requires much more 
computation effort than linear extrapolation. However, it 
has the advantage that if the time series has a slowly time 
varying structure, the AR equation can model it and provide 


accurate forecasts. 
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3. SSF and Conventional Feedback and Feedforward Control 


3.1 Introduction 

The conventional feedback and feedforward COlntoor 
systems are the oldest and still the most common control 
systems used in chemical process control. Their structures 
are very simple and the theories behind them are well 
understood. The purpose of this chapter is not to _ propose 
any new theory about these conventional control systems, but 
rather to identify some of the concepts used by them. These 
concepts will be shown to be intrinsic in more advanced 
control systems such as the Internal Model Control discussed 
in later chapters. The presentation is based on an 
intuitive, practical approach such as would be obtained by 
experience with these control systems. Another purpose of 
this chapter is to explore the relationship of SSF with 
these traditional control systems. It will be shown that 
the proportional-derivative feedback controller can be 
interpreted aS a proportional controller with a linear 
forecaster. In addition, the feedforward control concept 
has been extended in a scheme which uses SSF _ for 


compensating unmeasured disturbances. 
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3.2 SSF and Conventional Feedback Control 


3.2.1 Description of the Feedback Control System 

The block diagram of a conventional PID feedback 
control systemjis shown in Figure 3.1. neg OL OGh@EsC la Giaam 
has been simplified by omitting the measurement device and 
the fanel controlyelement since theiridynamics) are typically 


negligible-and/or incorporated in the process model. 


Ys ; 
SETPOINT + 


CONTROLLER 


af 
PROCESS OUTPUT 


Biguresssral Block Diagram of a Conventional Feedback 
Control System 


The principle of feedback control is the feedback of 
The CuUcpUL Variableytoecalculatethescontrol enror,se) which 


forms the input to an “error driven" controller. 
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Gh=SY. a Uy (37) 


The control error can arise from two sources: setpoint 
changes and disturbances. Servo control refers to 
eliminating the control error arising from setpoint changes. 
Regulatory control refers to eliminating the control error 
arising from disturbances. Since the error-driven feedback 
controller does not distinguish the source that gives rise 
to the control error, it cannot distinguish between servo 


and regulatory control. 


CONTROLLER MODES 

The conventional PID feedback controller has three 
modes of eontrol action: “proportional, integral), and 
derivative. Mathematically, the controller output in the 


continuous time domain is given by: 


1 d 
u(t) ==Koele(t)- +.- * fet) dtbeterg—e (t) ] (352) 
Tay dt 


The Simplest feedback controller has only the 
proportional mode. The problem with it is that it leaves 
offset in the output. The integral mode is added primarily 
to eliminate this offset. However, it also makes’ the 
response more oscillatory. The derivative mode is added to 


speed up the response of the system. 
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Equation (3.2) represents the theoretical PID feedback 
controllers. In practice, there are many variations. One 
variation which is of interest to this study is where the 
derivative mode acts on the process output instead of the 
CONtLOl error, PNeeDGactice, mihi cman: s a dOnemtOUmavyOl1dema 
"derivative kick" which is a sudden change in the controller 
output when there is a sudden change in the setpoint. The 


equation is 


. 1 d 
Wc me=eKcm te GC) it a i'eGl) dtetanrg—yt) | (sie) 
T j dt 


CLOSED LOOP TRANSFER FUNCTION 
The closed loop transfer function of the system in 


Figure 3.1 is given by: 


Gc Gp Gd 


ee ———— 
1gteGesGp int GcaGp 


The setpoint y, and disturbance & are the two inputs to 
the closed-loop system. The transfer function between y and 
y; determines the dynamics of servo control while the 
transfer function between y and £ determines the dynamics of 
regulatory control. Since the two transfer functions are 


different, the dynamics in servo and regulatory control are 
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different. Therefore, the use of a single controller in the 
feedback control syStem may not guarantee good performance 
in both cases. There may be advantage in separating the two 


controllers. 


CONCEPT 3.1 Servo control and regulatory control 
generally have different closed loop dynamics, it 
would be better to have separate servo and 


regulatory controllers. 


From equation (3.4), perfect control can be defined as 


requiring 


Gc Gp 
——__—_—_—_—_ = 1 (ols SYED) 
tat+aGe 7Gp 
and 
Gd 
———————- = 0 (390) 
i=teGcacp 


This is possible when the gain of the controller is 
very high or infinite. Hence, perfect feedback control can 
be obtained with an infinite gain controller. 

However, the gain of a feedback controller is usually 
limiteds byfostabilityseconstraints. According to the Bode 
criterion, a stable loop requires that |GcGp| < 1 at the 
crossover frequency i.e. when the phase angle is -180°. 


Therefore, the controller gain iS constrained and perfect 
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control in the sense of equation (3.5) is not always 


possible. 


CONCEPT 3.2 Perfect feedback servo and regulatory 
control can be approached by using a high gain 
controller provided that the high gain is within 


the stability constraints. 


3.2.2 Interpreting Proportional-Derivative Control as SSF 

One characteristic of feedback control is that it is 
remedial i.e. it corrects for error which is already present 
in the output. For systems which contain time delay or 
measurement lag, feedback control is not very good because 
the control action always lags behind the disturbance and 
tends to overcompensate. AWISeCISseretes dig@rralm control 
Systems have at least one period of inherent time delay due 
to sampling. The approach used to control time-delayed 
system is to have some form of prediction. In this section, 
it will be shown that the proportional-derivative controller 
contains a form of forecasting. 

Consider the structure of a proportional-derivative 
controller in a discrete system (Figure Bhd Fe The 


controller can be considered as composed of two parts: a 
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forecaster and a controller. The forecaster is given by 


(1 + ry —) e(t) (30) 


FORECASTER CONTROLLER 


Rigusegs. 2 Structure of a Proportional-Derivative 
Controller, 
For the special case where ry, = n*xAt, the expression 


forecasts e(k) exactly n sampling intervals ahead using the 


slope of e(k). 
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The controller 1s a proportional) controller) “with ‘gain 
Ke which acts on the forecasted error é(k+n) and equation 
(3.7) is equivalent to the equation of a linear forecaster 
as given by equation (25208 Therefore, the 
proportional-derivative controller can be interpreted as 
containing a linear forecaster of the contro] error. A 
one-step ahead forecaster is needed to compensate for _ the 
inherent sampling delay in discrete control system (Jacobson 
(1970) and Moore (1969)). Forecasting ahead for more than 


one step can be used to compensate for delays and/or process 


lags. 


CONCEPT 3.3 The proportional-derivative controller 
in Figure 3.2 can be interpreted as a proportional 
controller with a linear forecaster on the control] 


error. 


Consider now the feedback controller described by 
equation (3.3) where the derivative mode acts on the process 
GuEpUL rather #thany the controlwerror. "A “block diagram oF 
the control system with a proportional-derivative controller 
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FEEDBACK CONTROLLER 


Figure 3.3 A Re-structuring of the PD Feedback 
Control System 


The controller in the feedback path is the same PD 
BeedbaGkemiCOnonOll Gr acoS, atic t met eels Olen se Shere The only 
di fherenceeahS suthat aches SinputeestoOu. the controller is the 
process output y. Therefore, following the same argument 
hore the sleedbeck controller insFigure 3.2, 1t.can. be shown 
that the feedback controller in Figure Sa can be 
interpreted as containing a linear forecaster on the output 


since 
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The controller gain Kc acts on the forecasted process 


CuUEpUL Ry (krie 


CONCEPT 3.4 The proport ional-derivat ive controller 
in Figure 3.3 can be interpreted as a proportional 
controller with linear forecasting of the process 


output. 


Another feature in Figure 3.3 is that there are two 
controllers in the system. The one in the "Setpoint path" 
can be regarded as the servo controller and the one in the 
"feedback path" can be regarded as the regulatory 
controller. As mentioned previously, the dynamics in servo 
control can be quite different from that in regulatory 
control. Therefore, with this separation, it is possible to 
design different servo and regulatory controllers hte 
justified by improved performance. The implementation of 
two controllers is easily done in computer control system 


since all that is required is software. 


3.3 SSF and Feedforward Control 


3.3.1 Description of Feedforward Control 
Feedforward control is used primarily to compensate for 


a disturbance before it upsets the process. The block 
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diagram of a conventional feedforward control system is 


Shown in Figure 3.4. 


PPe-b aBeCONTROLEER PROCESS 


Figure 3.4 Block Diagram of a Feedforward Control 
System 


The disturbance that affects the process is measured as 
it enters the process. Based on this measurement, the 
feedforward controller generates appropriate control action 
to compensate for the effect of the disturbance so that’ the 
output 1s not affected. Note that a key feature of 
feedforward control is derived from the fact that the 
disturbance is known and corrective action is taken before 
its effect propagates through the opens to affect the 
output, This time advantage gives the controller time to 


take action to compensate for the effect of the disturbance. 
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CONCEPT 3.5 Feedforward control makes use of the 
time advantage from knowing the disturbance and 
taking corrective action before its effect 


propagates through the process. 


The objective of the feedforward controller can be 


described by the following equation. 


£) GaesF SG xeGp = (0 (3729) 


From equation (3.9), the structure of the perfect 


feedforward controller is obviously 


Gx =-—— (53 a0) 


The controller works on the principle of cancellation 
of transfer functions. For perfect control, i.e. exact 
cancellation, the disturbance and process transfer functions 


must be known exactly. 


CONCEPT 3.6 In conventional feedforward control, 
the perfect controller transfer function is a 
product of the load transfer function and the 


inverse of the process transfer function. 
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The combination of Gd and Gp to form Gx often results 
in simplification in the structure of the feedforward 
controllers For example, consider a first order process 


with a first order disturbance tranfer- function: 


Kp 
Gp = ——— Ossie) 
%, £5: 4+ ol 
Ka 
Gd = ——— (3b) 
TiS ee | 
Lint = ts enen 
Kd 
Gx = - —— = - Kx 3 ea12) 
Kp 


The feedforward controller is Simplified ZO a 
proportional controller. 

Note that the perfect feedforward controller requires 
the inversion of the process transfer function. This could 
be a problem if the process transfer function contains’ some 
non-invertible parts such as a time delay or unstable zeros. 
One solution is to find an approximate inverse for the 
process which neglects these non-invertible parts. In 
practice, most feedforward controllers are proportional 
controllers Or have simple "lead-lag" type dynamic 


compensation. 
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3.3.2 Feedforward Control Using SSF 

The disadvantage of conventional feedforward control is 
that it is designed only to compensate for measured 
disturbances. However, some of its principles can be 
applied to compensate for unmeasured disturbances ae 
estimates of these disturbances can be obtained. It will be 
shown in the following how SSF can be used as part of a 
scheme to compensate for unmeaSured disturbances in a 
feedforward manner. 

Figure 3.5 shows a block diagram of the control system 
to be described. The first step in the control scheme is to 
obtain an estimate of the effect of the unknown disturbance 
on the current output. This is achieved by introducing a 
process model which estimates the process output at every 
Sampling instant. Not much will be said about the process 
model at this point except that its function is, given 
values of the process input variables, to estimate the 
process output as accurately as possible. The type of model 
used is not of concern here. 

The difference between the actual output and the 


estimated output is the estimation error or residual r(k). 


ECk) =a yk) eeomy Ck) (ints 


If the model is exact, the residual is equal to the 


effect of the unknown disturbances y,(k). If the model is 


not exact, the estimation error will contain some modelling 
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error as well as the disturbance. Note that in the ideal 
Casereithe @residual Pisthan pestimate softi the unmeasured 


disturbances. 


CONCEPT 3.7 A residual can be generated’ by 
subtracting the predicted output from the measured 
output. This residual 1S) "ean estimate of 


disturbances and/or modelling error. 


In order not to complicate the discussion for the time 
being, it will be assumed that the model is exact and 
mek), =iy, (kde 

It is not sufficient just to know the effect of the 
unknown disturbances in the current output. For feedforward 
controls the time advantage in knowing the future 
disturbances (CONCEPT 3.5) must be provided. This is 
achieved by introducing a forecaster which forecasts the 
future residual r(k+A) based on the past values r(k), 
r(k-1), ... . The type of forecasting technique used is also 
not of concern here. It may be any one of those mentioned 
in Chapter ~lwor The main criteria for the choice of the 


forecaster are accuracy and computational effort. 


CONCEPT 3.8 The time advantage required for 
feedforward control can be generated by 


forecasting the estimated disturbances. 
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The controller design is analogous to the feedforward 
controller described in the previous’ section. Given the 
estimated disturbances, the controller that will generate 
the exact compensatory action igs a predictive controller. 


Its form is 


Co (2514) 


Since thetinpupetosethistcentroller is ty esrather ithan?é?} 
Gd is not present in equation (3.14). The minus sign which 
was in equation. G3 £10) Stbute® noti iimnerc3 914) Behas @ibéen 
incorporated in the summing junction in Figure 3.5. 

This controller also works on the principle of 
cancellation of transfer functions. Fone periectescontroly 
the process transfer function must be known exactly i.e. 


Gm = Gp. 


CONCEPT 3.9 Perfect feedforward control based on 
estimated disturbances requires the use of an 
exact process model for prediction of the process 


output and an exact inverse for the controller. 


The controller design determines how far ahead the 
forecaster has to estimate. This is actually related to the 
time delay of the process. Assume that the process contains 


a time delay of A sampling periods including the unit delay 


¢tt £&) 


ni 

2 vs net Sas 4g ei Sansa gh ws ‘suqani off eonte: 
faidw. apie: suite “ai?! eeae Reape wi tsgesoqg zonial 
aged odd COTTE} ont Sgon tud (07,2) nottsups (ae ” 3B _ 

ace” abuse nal nekasaay wrt ae ni 7109 

a vabersdtva ale 8 ekzow ae ls , Hises rete” aici 
evans rs@kaey tot agian: teRebe= 2 20 neisall 

ek MADRE RO) | sie if ele 9 oon —ee phienass a 

! 


a 


ate ADR 2 2.4; 


20 


due to sampling. When the process transfer function is 
inverted to form the controller, the delay becomes a 
prediction. Therefore, the forecaster has to estimate A 
sampling periods ahead to provide the disturbance estimate 
for the controller. 

In the inversion of the process transfer function, a 
problem is encountered if the process transfer function 
contains unstable zeros. The controller will be unstable if 
the inversion is done. Therefore, the above controller 
design is restricted to stable-inverse processes only. For 
process with unstable zeros, an approximate inverse may be 
used by inverting only the stable and invertible part of the 
process. This problem will be addressed again in Chapter 
Four in the discussion of IMC. 

A simulation example will be given in the following as 
an illustration for the system in Figure 3.5. Consider a 


first order process described by the discrete ARMA equation. 
yee OT OvCkK lect 0). 60 (Kos) meter (Ket Sf Ks @.1 99) 


The effect of the disturbance on the output is y,(k) 


and is assumed to be 
vVe(kjlevtdk wimOeSét kone (3216) 


Assuming the disturbance is unknown, the best 1/0 


process model is 


+ 28. 
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y(k) = 0.8y(k-1) + 0.6u(k-1) (37917) 


Two forms of one-step ahead forecasters are considered 
here. The first one is a linear forecaster described by the 


following equation 


mk type—eer CK eter kK) Ser (k= 1)) (34.8) 


The second one is an adaptive autoregressive forecaster 


of the form 


HOGkK+ br =Gage (ky) € fa gr(k—-1 Pisuw er CkS2) (34199 


The parameters a;, a@2 and a3 are identified on-line by 
a parameter estimation scheme. 

The same process model (equation (3.17)) has been used 
to design both Gc and Gx as per Figure 3.5. The controller 
output is derived from equation (3.7) by replacing y(k) with 


Vek ie oa CK +1 )e, 


u(k) = Ly. (k+1) - £(k+1) - O-8y(k)] / 0.6 (3*209 


In?.the Simulatronsec(FiguressacGatoisel 1)ipethetprocess 
was initially at steady state and an arbitrary disturbance 
Sequence was introduced. The process was first simulated 
under feedback control using well-tuned controllers in 


Figure 3.6 andees. ss In Figure 3.6, a proportional 
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proportional-derivative controller of the type shown in 
Figure oie was used and an improvement in control 
performance was obtained. This improvement can be 
attributed to the inherent forecasting of the 
proportional-derivative controller which tends to offset the 
process lag. 

EPNLExrgquressecs tohe3.11® thet@process was underetthe 
control of the proposed feedforward control system in Figure 
3.5. The process model was perfect as given by equation 
(3.17) while the controller was an inverse model controller 
as given by equation (3.20). Figure 3.8 shows the system 
response for the case of perfect forecasting of the 
disturbance effect, the process output is not affected at 
alle This is the best and ideal case. No forecaster was 
used in Figure 3.9 and there waS no compensation for the 
disturbances. This is the worst case for the system which 
is under inverse model control. 

Figure 3.10 and 3.11 show the response of the system 
under more realistic conditions. The linear forecaster of 
equation (3.18) was used in Figure 3.10. An adaptive 
autoregressive forecaster (equation (3.19)) was used in 
Figure 3.11. The parameters in the adaptive forecaster were 
identified by the projection algorithm used in the Adaptive 
Predictive Control System. In this simple example, the 
linear and adaptive forecasters appear to give similar 


control performance which is better than that given by the 
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proportional feedback controller in Figure 3.6. However, 
the input sequence for the adaptive forecaster is smoother 
than that of the linear one. 

The response of the proportional-derivative controller 
(Figure 3.7) is different from that of the SSF system with a 
linear forecaster. The reason is that the 
proportional-derivative controller useS a proportional 
controller whereas the SSF system uses a predictive 
eontrollen. 

This simulation example shows that the use of SSF can 
improve control performance with respect to disturbance 


compensation. 
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Figure 3.6 Disturbance Compensation by a Proportional 
Controller (Kce= 17) 
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Figure 3.7 Disturbance Compensation by a Proportional- 
Derivative Controller (Kc=1.7, Ky=1.0) 
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Figure 3.10 Disturbance Compensation by a Predictive 
Controller with Linear SSF 
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Figure 3.11 Disturbance Compensation by a Predictive 
Controller with AR SSF 
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3.3.3 Comparison of Conventional FF and SSF FF Control 

SSF has been proposed for the compensation of 
unmeasured disturbances. It is a modification of the idea 
of conventional feedforward control. The differences 
between conventional FF control and SSF FF control are 


summarized in the following table. 


Conventional FF FF using SSF 


1. compensates for 1. compensates for estimated 
measured disturbances residuals, e.g. unmeasured 
only disturbances and modelling 

errors 

2. measures disturbances 2. estimates disturbances 
directly as they based on measured process 
enter process output and estimated output 


generated from a process 


model 
3. has an inherent time 3. needs to extrapolate 
advantage by measuring SOpCONtLOWNaACt1 ON Uk) 
disturbances as they can compensate for y,(k) 
enter the process 
4, perfect controller is 4, perfect controller is 
Gx = - Gd / Gp Gx =m Gp 


(note minus sign in summing 


junction in Figure 3.5) 
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The implementation of conventional FF control requires 
hardware for the meaSurement of the disturbances. On the 
other hand, the implementation of SSF FF control ina 
computer control system requires computer software only. 
Therefore, any computer control system can add feedforward 
control using SSF instead of measurement hardware. However, 
the performance of SSF FF depends on the quality of the 
process model and the forecast. The availability of a good 
process model and @eforecaster® is, thereforep? val’ major 


consideration in using SSF FF. 


3.4 Conclusions 

CONVENTIONAL FEEDBACK CONTROL 

1) Proportional-derivative control can be interpreted as 
proportional control with a linear forecaster. 

2) Forecasting of the output signal can improve control 
performance when a system contains time delay due to 
sampling and/or process lags. 

3) Perfect control in feedback systems can, in theory, be 
achieved by an infinte:gain controller. However, this 
is impractical because the gain is limited by contraints 
imposed by stability requirements and physical 
realizability. 

4) Separation of servo and regulatory control functions may 
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FEEDFORWARD CONTROL BASED ON ESTIMATED DISTURBANCES 


1) 


2) 


3) 


A control structure which consists of a process model, a 
forecaster and an inverse model controller has been 
shown to be capable of compensating unmeasured 
disturbances in a feedforward manner. Simulation 
results show improvement in performance of the SSF 
System over conventional feedback control. 

Perfect feedforward control is achieved by the use of an 
exact model inverse controller. The principle behind it 
is the cancellation of process transfer function i.e. 
Gc Gp = 1 and hence Gc = Gp"'. 

Feedforward control using SSF can be implemented as part 
OfACtheascontrolepcaliculatronceiny anyhacomputesicontrol 
System. It does not require the addition of hardware to 
measure the disturbances. All that is required for 


implementation is computer software. 
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4, SSF and Internal Model Control (IMC) 


4,1 Introduction 

Internal Model Control (IMC) was proposed by Garcia and 
Morari (1982) as an attempt to unify several different 
control schemes (Smith Predictor, Inferential Control, Model 
RiIgoumehmicws Control,» Dynamic MatrurxesCcontrol, ete.) 2.5) The 
interest in formulating IMC is partly aroused by the success 
of the two control schemes developed in industry: Model 
Algorithmic Control (MAC) and Dynamic Matrix Control (DMC). 
They were developed on a heuristic basis and have shown 
excellent performance in industrial applications. 

MAC was developed in France by Richalet et al (1978) 
and its theory was extended by Mehra et al (1980). DMC was 
presented by Cutler and Ramaker (1980) of Shell Oil in the 
United States. A constrained multivariable version of DMC 
was applvedietom ea’? catalytic® eracking Punite byePrecteand 
Gillete (1980). 

The other work which stimulated the proposal of IMC is 
the inferential control structure proposed by Brosilow 
(1979). Ind a SStudy#%of thevclassical’ Smith’ Predictor and 
inferential control, Brosilow proposed the Inferential Smith 
Predictor structure together with a set of controller design 
procedures. The IMC structure can be considered as an 
extension of this Inferential Smith Predictor structure 


which can also be shown to include MAC and DMC as particular 


cases. 


42 


ae 2 | itseyd 
neSspoue os ee baweily ictal 3 Wer prbselemro3 esusaai 
halam + ygasecbnt n> bem. eabe gahtnte Dodie ~ ows, sit 30 
ose) fortaed Sioa aacci S:alABin towsn0d ) | 
wmoda sved bus efkad sitziiced Ss nm Beqolaved: — 
atm igesligds’ tai? 2onme nt Sonamyoizeq 4 T 
(eter) te ts selatsit wa ae jbsgetaver: now at 7 
See Mer < GREP? Steer bata oS pebdesee bow 
ats di Dio Liane Ye Cane) en meee 146 
War ib pnelteiv sided pe eale}iereo: & | 
fis sec? Yo phew Pie heer 31 fdies B oF boteeash= 
ae | + (O8tY 9: 
nt OM? to Ianogetg eff Beter umise ashdw S20~ teto™ edt - 
wolieexd yc beaoqouy suc puise idzines ae 
fas Soto lSaxd ‘dadme bestersts- shgp-do- youte 6 aD 
cate, sail adiiel ial Aseoqoxg wiles slossaer 


ab 


ees a 


_ 


43 


The purpose of this chapter is two fold. The first one 
is to interpret the IMC structure according to the concepts 
of feedback and feedforward control as summarized in Chapter 
Three. The second one is to investigate the incorporation 
OPO SHae in CiNewmiMGCesStructiine. Simulation examples are 
presented to illustrate the properties of the modified 


system. 


4.2 Description of Internal Model Control (IMC) 
BASIC IMC STRUCTURE 

According to Garcia and Morari (1982), the IMC 
structure can be described in three levels of complexity. 
The first level is the basic IMC structure which is shown in 
Figure 4.1. 

The basic IMC structure is obtained from a conventional 
feedback control structure by adding and subtracting an 
CUEDUL, EY mpe aSashown in Figures4 ie ltsisvequivalentitora 
conventional feedback structure since the two paths 
incorporating Gm Cancede= Outs Theoretically, the IMC 
structure does not impose any reStriction on the type of 
process model used. In practice, an impulse response model 


is often used. 
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Figure 4,1 Basic IMC Structure 


bom Faqurese4ciue GC GZ) ei SmethnesIMGycontroilerts sitans 
related to.the conventional feedback controller C(z) by the 


following equation. 


Giz) 
Gc(z) = ———————___——_ (aah) 
1 + Gm(z)C(z) 


The claimed advantage of the IMC structure over the 
conventional feedback control structure is that Gc(z) is 
Sasierme to sGeSlgne st hate G( zi). This is because the design 
procedure for Gc(z) can be easily established from the IMC 


structure. Moreover, robustness can be included as a design 


- 
a 


45 


objective ina very explicit manner. This is partly due to 


the special form of the feedback signal d(z). 


a(z)m= [SeeueCGp( zdriGm (a) Gd tz) elie’ dbz) (4.2) 


If the model is perfect, the feedback signal d(z) is 
equal to the disturbance d(z). If the model is not perfect, 
d(z) will contain some information about the model-process 
mismatch. By modifying d(z) appropriately, robustness can 
be obtained. 

The following transfer functions can be obtained for 


the structure in Figure 4.1. 


Ge (z) 
u(z) = —————___ (y,, (z) -d(z) ) (4.3) 
1 + Ge(z)[Gp(z)-Gm(z) ] 
Gp(z)Gc(z) 
y(z) = d(z)+————___ (y.. (z) -d(z)) (4.4) 


1 + Ge(z)[Gp(z)-Gm(z) ] 


Equation (4.4) indicates that perfect control requires 


1 
Gce(z) = ———— (4.5) 
Gm(z) 


Because Gc(z) has to be stable, Gm(z) is factorized 


into two parts. 


Gm(z) = Gm,(z) Gm_(z) (4°56) 
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Gm,(z) contains the non-invertible part (time delay and 
unstable zeros) and Gm_(z) contains the invertible part. 
According to Garcia and Morari (1982), from optimal control 
theory and assuming Gm(z) = Gp(z) and Gp(z) is stable, the 
controller that minimizes the sum of the square of control 


errors is 


Gce(z) = ————— (4.7) 


COMPLETE IMC STRUCTURE 

It has been shown in equation (4.5) that the IMC 
controller is an inverse model controller. If the inverse 
is exact, the controller gives perfect control. However, 
very often the exact inverse cannot be obtained and an 
approximate inverse has to be used. This brings up the need 
LOG additional features for the IMC. structure which 
corresponds to the second level of complexity in the IMC 


Structure. Figure 4.2 shows the complete IMC structure. 
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Figure 4.2 Complete IMC Structure 


A filter has been added in the feedback path. Its 
Function as = tom account |fory the, model-processm mismatch 
reflected in the feedback signal d(z). Tpeecan also 
compensate for certain types of disturbance dynamics. The 
characteristic equation with reference to equation (4.4) is 


changed to 
Ge-'(z) + F(z) [Gp(z)-Gm(z)] = 0 (4.8) 


For a given model-process mismatch, F(z) can be 
designed to ensure stability i.e. equation (4.8) has stable 


roots. 
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Another feature in the complete IMC structure is the 
reference model" R(z) for “the ‘setpoint. It provides a 
reference trajectory which has certain desired dynamics. It 
can also desensitize the process with respect to modellina 
inaccuracies by "Shaping" the desired value, ya, seen by the 
Coneroller, Gc. 

In summary, the IMC controller design procedure for the 
ideal case involves selecting an inverse controller which 
gives perfect control. However, if an inverse controller is 
not feasible, an approximate inverse model controller is 
used anda filter is added to improve robustness. Finally, 
a reference model can be added to reduce sensitivity to 


modelling inaccuracies. 


IMC WITH PREDICTIVE CONTROLLER 

The third level of complexity in the IMC structure 
involves the design of a predictive controller. This is one 
of the ways of obtaining an approximate inverse of the 
model. This approach has been found to possess many 
desirable properties. The formulation of the predictive 
controller is given as follows. 

The predictive control strategy considers a desired 
output trajectory yu(k) over a horizon of P sampling times 
into the future. Then the sequence of control actions u(k), 
UK? ee oie UKeP op emewhel Cl Kin om HeRmGUCTen Db amcamo ling 
instant, is calculated so that the predicted output 
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closely as possible. The prediction is calculated by the 
Hodes Oe Dem—E I DUtS Up CORK tla mpmwuel Caged Smtne pune 
time delay, and output up to k. If the process model, which 
is used for the prediction and the control calculation, is 
not exact, the process output will deviate from the desired 
trajectory. Therefore, it is preferable to implement only 
the present input u(k) and to resolve the problem again at 
k+1 with -the measured output y(k+1) as the new starting 
point. Figure 4.3 illustrates one step in this moving 


horizon problem. 


Figure 4.3 One Step in the Predictive Control Problem 
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At each sampling instant k, the following 


problem is solved. 


P 
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subject to 


G(k+r|k) = ym(k+r) + A(k+r|k) 
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where P is the horizon (P 2 1) 


ya(kt+r+i) is the desired trajectory 
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general 


Goatees 2) eectmcr, oem ten CCE) EC ak ope) 


y,? are time varying weights on the output error 


B,;? are time varying weights on the input 
M is the input suppression parameter which 


Specifies the number of intervals into the 


future during which u(k) is allowed to vary 


y(k+r|k) is the predicted output 
ym(k+7) is the output of the internal model 


d(k+r|k) is the predicted disturbance. 


The simplest prediction of the disturbance is to set it 


equal to the residual at the present time. 
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A(k+7r|k) = A(k) formal lbkwe.dk (4.10a) 
A(k) = y(k) - ym(k) (4. 10b) 
M, Pros, tro thiareecthem tunangMiparameters: gofl the 


algorithm. They have a direct influence on stability and 
dynamic response. 

The complete IMC structure with predictive controller 
is shown in Figure 4.4. Details on the computation of the 


control law were presented by Garcia and Morari (1982). 


ya(z)} OFFSET SETPOINT 
COMPENSATOR PREDICTOR 


Figure 4.4 IMC Structure with Predictive Controller 
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4.3 Discussion of IMC 

Garcia and Morari (1982) showed that the IMC. structure 
is related Low Optimal’ = Control, “the “Smith =Predictor, 
Inferential Control, Model Algorithmic Control and Dynamic 
Matrix Control. In this section, the IMC structure will be 
interpreted according to the concepts of feedback and 
feedforward control as summarized in Chapter Three. The 
issues of accuracy and invertibility of the process model 
will also be explored by examining the filter, reference 


model and predictive controller in the IMC structure. 


4.3.1 Feedback Interpretation of IMC 

The first step in the IMC controller design is to 
select a controller which gives perfect control. It has 
been determined from equation (4.4) that the perfect IMC 
controller is the inverse of the process model. In order to 
understand how the inverse model controller achieves perfect 
control, the equivalence of the relationship between the 
feedback control system and the IMC structure will be 
examined. 

Equation (4.1) gives the relationship between the 
feedback controller C(z) and the IMC controller Gc(z). The 


equation can be re-arranged into the following form. 


Ge (z) 
CARs ia?) 


C(z) = 
1 - Gm(z)Gc(z) 
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Substituting ’ the perfect IMC controller, 
Gc(z) = Gm-'(z), into the equation, the denominator of the 
expression becomes Zero and C(z) becomes infinite. 


Therefore, the IMC inverse model controller is equivalent to 
an infinite gain feedback controller. Perfect control in 
IMC, as in feedback control, is achieved on the basis of 
infinite gain (Concept 3.2). 

An alternate way of obtaining the same result is 
through the block diagram. Figure 4.5 shows an equivalent 


IMC structure which is drawn differently from Figure 4.1. 


Figure 4.5 AnaEQuivalent.Porm of the siMC Control ler 


The IMC structure of Figure 4.5 becomes a conventional 


feedback control structure if the positive feedback loop is 
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included in the controller block. Using the perfect IMC 
controller, Gc = Gm~', the positive feedback loop becomes an 
infinite gain controller in a feedback control structure. 
therefore, (Piromi@aigstructural \eporntpeolfercviews, IMC is 
equivalent to feedback control and perfect control is 
achieved through the principle of infinite gain. 

In a feedback system, the use of high gain is limited 
by stability and performance considerations and perfect 
control may not be attained under such restrictions. 
Similarly, the IMC inverse model controller can only attain 
perfect control when the model is perfect and the inverse 
model controller is realizable. The reason for this is as 
follows. 

For discrete systems, there is at least one sampling 
period of time delay. Therefore, the IMC inverse model 
controller will contain at least one period of prediction. 
This means that the controller input, ya - dG, must be known 
at least one period of time ahead. The setpoint yg can _ be 
specified ahead Offetimes However, if the model is 
imperfect, the feedback signal @ will contain model-process 
mismatch and disturbances. Future values of da cannot be 
known since the model-process mismatch is dependent on _ the 
current process input which needs to be calculated using 
future values of d. Therefore, perfect controm is 
impossible with an imperfect model. ingvaddition;, 
realizability of the perfect controller requires that the 


process be minimum phase. 
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Since future values of d are required in the control 
calculations, they are approximated by the current value of 
Gd in the IMC structure. This approximation represents a 
deviation from the ideal condition of perfect control. The 
extent of this deviation depends on the model-process 
mismatch contribution to the feedback Signal. The 
contribution of the mismatch decreases as the process model 
becomes more accurate. When the process model is perfect, 
the feedback signal d contains disturbances only. For the 
case of no disturbance, the controller input is equal to the 
setpoint which is known exactly and perfect servo control 
can be obtained. When disturbances are present, perfect 
regulatory control can only be obtained if the disturbances 
can be forecasted accurately. Hence the interest in 


combining SSF with IMC. 


4.3.2 Feedforward Interpretation of IMC 

An alternate interpretation of IMC is based on the 
feedforward control structure using SSF (Figure 3.5). This 
interpretation is appropriate for explaining the need for a 
pertect! process #imodeli® andsa mone Sraccurately gmeidiectsim the 
design philosophy of IMC. 

The IMC structure (Figure 4.1) is similar to the 
estimated feedforward control structure (Figure 3.5). Both 
of them make use of a process model to generate a _ residual 
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control, the controller has to be equal to the inverse of a 
perfect process model. It has been shown that perfect 
control in feedforward control is achieved through 
cancellation of transfer functions (Concept 3.9). Because 
of the similarity in structure between IMC and feedforward 
control, perfect control in IMC can also be interpreted as 
being achieved through cancellation. This @trs? ‘avcdimore 
appropriate interpretation than the feedback interpretation 
when the process model is perfect since the feedback signal 
is equal to the disturbance and the system is effectively 
open loop. It should be noted that the forecaster present 
in the feedforward control structure is not obvious in the 
IMC structure. However, some form of prediction is implicit 
in the inverse model controller. 

The IMC structure can be interpreted as being designed 
for open loop feedforward control using the inverse of a 
perfect process model as the controller and an estimated 
disturbance d(z) (which is equal to d(z) when the process 
model is perfect). A perfect model removes the problem of 
instability arising from the feedback control structure. 
The perfect process model is also used to generate the 
residual signal which is used for regulatory control in a 
feedforward manner. This is the ideal IMC design 
philosophy. When the model is imperfect and/or 
non-invertible, then feedback is present, instability could 
result and other features such as the filter, reference 


model and predictive controller are necessary. 
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4.3.3 Filter and Reference Model 

The filter and the reference model can be considered as 
part of the controller in IMC. When the filter and the 
reference model have the same transfer function, they can be 
combined with the controller to form an equivalent system 
with only one control block. The advantage of separating 
part of the controller out to form the reference model in 
the setpoint path and the filter in the feedback path is to 
provide a means to have different servo and regulatory 
controllers. The reason for it has been discussed in 
Chapter Three (Concept 3.1). 

The function of the reference model is to generate a 
reference trajectory from a given setpoint for the system to 
track. A reference trajectory can often be made more 
appropriate than an external setpoint because it usually 
requires less extreme control effort to track and it may 
possess certain optimal characteristics. In addition, it 
can desensitize a system i.e. make the system respond 
slower. 

The filter, on the other hand, can be used in a number 
of ways. The first use of the filter proposed in the IMC is 
to compensate for the model-process mismatch. From equation 
(4.8), it can be seen that the filter can be designed to 
ensure stability and shape the dynamic response of the 
system for a given model process mismatch. Brosilow (1979) 
has discussed the use of a filter to compensate for 


modelling errors in time delay, process gain and process 
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time constant using the example of a first order process. 
He@ showed Dthat@tarstirst torder lag filter with an 
appropriately chosen time constant will provide satisfactory 
compensation in most cases. 

The filter can also be used to compensate for certain 
classes of disturbances. A typical example is that of a 
low-pass filter which is used to filter out high frequency 
disturbances or noise. This use iS very important in 
providing smooth measurement data to the predictive 
controller. 

The filter can also be interpreted as generating a 
reference trajectory for the disturbance signal in the 
feedback path for regulatory control. This helps to reduce 
extreme excursions in the process input. It slows down the 


response of the system and improves robustness. 


4.3.4 IMC Predictive Controller 

The perfect IMC controller is an inverse model 
controller. “Therefore, the invertibility of the “model is 
very» important “in’the controllersdesign. For models which 
contain non-invertible parts, an approximate inverse has to 
be used. There are two ways of obtaining an approximate 
inverse of the model. The Gfarstthonereisettoe® invertes the 
invertible part and approximate the non-invertible part. 
The second one is to approximate the inverse of the whole 


model. The predictive controller in IMC inverts the process 
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model using the second method. With appropriate choice of 
the tuning parameters, the inverse can be made exact or 
approximate. The objective in the controller design is to 
find a stable inverse of the model which gives good control 
performance. The physical significance of the tuning 
parameters and their effects on the controller 
characteristics have been discussed by Garcia and Morari 
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4.4 SSF and IMC 


4.4.1 Description and Interpretation 

The IMC structure requires the use of an inverse model 
controller. In discrete control systems, an inverse model 
controller implies that some form of prediction is inherent. 
In the IMC structure, this prediction requirement means that 
future values of the controller input, 1.e. Setpoint and 
feedback signal, are required for the control calculations. 
To simplify the discussion, the model is assumed to be 
perfectieso thaterthe «feedbackstsrgnal ©isesequaliyto the 
disturbance. 

Consider the IMC structure with predictive controller 
in Figure 4.4. The need for future setpoints does not 
present any problem since they are simply assumed to be 


specified by the operator. However, values of future 
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disturbances have to be generated by some means. The 
predictive controller proposed by Garcia and Morari (1982) 
uses the present disturbance in the place of the future 
disturbance (Equation (4.10)). There is no forecasting done 
to obtain the future disturbances. 

From another point of view, this can be interpreted as 
being equivalent to the first forecasting method described 
in Section 2.3. It uses the current value as the forecast. 
This is justified if the disturbance is stochastic and 
therefore difficult to forecast. However, in cases where 
the disturbance has some structure, it is often possible to 
forecast future values with a reasonable degree of accuracy 
by using information from the past values. This is exactly 
the idea of the SSF. It fits in the IMC structure and is 
used to forecast future disturbance values based on the past 
values. Figure 4.6 shows a schematic diagram of an IMC 
Structure with SSF. 

The forecaster can be interpreted as part of the 
inverse model controller in this context. It is added to 
make the inverse model closer to the exact inverse. The 
inverse model controller is then made up of three blocks: 
Gontcoller se stulLeiweands forecaster. Conceptually, the 
functions of the three blocks can be separated. The control 
block contains the inverse of the invertible part of the 
model. The forecaster contains the inverse of the time 
delay of the model. The filter contains a reference model 


which generates a reference trajectorvyetor better controle 
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Another interpretation of the forecaster is that it is 
a model of the unknown disturbances since it generates 
foOrecastsGof future disturbances. This > is “particularly 
appropriate if an adaptive AR forecaster is used. The 
accuracy of the forecasts depends on the forecaster 
structure and parameters. When there is doubt about the 
forecasted disturbance, the controller may want to 


compensate for only part of the disturbance. This idea is 


Similar to the concept of "relaxation factor" in numerical 
methods. inththesawcase, Thhedi filter, tewhich? contains’ @a 
reference model, can function as the "relaxation factor". 


By proper choice of the filter, good control can be obtained 
even if the forecaster does not forecast very accurately. 
Intuitively, the forecaster will forecast well when the 
disturbance is well structured and has_ slow dynamics 
relative to those of the process. It will forecast poorly 
when the disturbance has fast dynamics and no structure. In 
view of this, another use of the filter block will be to 
filter out the fast and random disturbances or noise. A low 
pass filter is well suited for this function. Then the 
forecaster can work on the slow filtered disturbances. In 
cases where the process model is a low order approximation 
of theshigh order process, the modelling errors which® are 
the fast modes of the process will be part of the feedback 
signal. Since the important process dynamics are generally 
Slow, wsthevlow passeftlter can alsomrilter out errorssduesto 


these fast modes. When the filter block is a low pass 
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filter, then for purposes of interpretation, the order of 
the forecaster and the filter should be reversed in 


Figure 4.6. 


4.4.2 Illustrative Examples 

A simulation example will be given to illustrate the 
use of SSF in IMC and to explore some of the properties of 
the system. A second order process without a process delay 
is used in the example. The discrete transfer function of 


the process is 


y(z) OeiS8ze' eh0eH0 1zer 
Ca a a ee GS be 
u(z) Trimet, 08 z6cnsn 02728 Szee 


The 8th order impulse response truncated model is 


y(k)e= 0.158u(k-1) + 0.264u(k-2) + 0.226u(k-3) 
+ 0.156u(k-4) + 0.096u(k-5) + 0.053u(k-6) 
be Or 02g u( kei \ereO> ONGutks 8) (4.13) 


The following set of parameters are used for the IMC 
predictive: controlder, inh ald@ethemicasesseandiethemeinpuc 


suppression parameter M is set individually in each case. 


baz Ni=l8 (4. 14a) 
y,Gety2 =) .2ho=eyen= Ul (4. 14b) 
B, = B2 = = BP, = 0 (4.14c) 
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Two forecasting techniques are studied. The linear 


forecaster as given by the equation 

Buk bbs rn Ok deta Gr Gk ed Gri Ck 1))s) (4.15) 
and the AR forecaster: 

mC ke )eeia inek) it age (k=t)et+ agrntk—2) (4.16) 


where a;, a@2 and a3; are identified on-line by the APCS 
projection algorithm (Equation (2.9)). For the simulations 
in Figures 4.7 to 4.10, the process is subject to a 
disturbance sequence which consists of four types of 
disturbances: ramp, exponential decay, sinusoid and step. A 
perfect process model is used for prediction and control 
calculations. The input Suppression parameter M for _ the 
controller is set to be equal to N. For this particular 
case where P=M=N, the number of unknowns is equal to the 
number of equations. An exact solution, rather than a least 
square solution, is obtained and the controller is an exact 
inverse of the process. Control action is calculated at 
every sampling instant. In cases where SSF 1s used, a 
one-step ahead forecast of the residual is calculated and 
used as the value for all the future residuals. The reason 
for this will be explained later. 

Figure 4.7 shows the process response under IMC control 


with no SSF. A constant offset is present during the ramp 
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disturbance. This indicates that IMC contains some form of 
integral control action. Figure 4.8 shows the response of 
perfect control as a result of using a perfect forecaster. 
Note that the forecast need only be perfect for the one-step 
ahead residual and not for the future ones. In the control 
calculation, an exact solution is obtained due to P=MEN. 
With control action being unconstrained, any setpoint change 
can be accomplished in one step. Since the calculation is 
done at every sampling instant, only the one-step ahead 
forecast of the residual is important. The future values 
has no effect on the implemented control action. This is 


why only a one-step ahead forecast is calculated. 
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The linear forecaster of equation (4.15) is used in 
Figure 4.9 and the AR forecaster of equation (4.16) is used 
in Figure 4.10. Comparing Figure 4.7 with Figures 4.9 and 
4.10 shows that the addition of SSF in IMC does improve 
regulatory control performance for smooth disturbances. SSF 
helps to eliminate the offset during the ramp disturbance 
and the non-linear disturbances. However, as expected, the 
forecasters do not help in the case of the step 
disturbances. Comparing the linear forecaster and the AR 
forecaster, the linear forecaster seems to be superior in 
performance for the smooth disturbance in this example. 
However, it gives rise to a more oscillatory input sequence 
and is very poor for step disturbances. 

Since noise iS an integral part of most measurement 
Signals, the following simulations will help to explore the 
effect of noise on the functionality of SSF. A noisy 
disturbance sequence is generated by adding random noise to 
the smooth disturbance sequence used previously. The 
conditions for the process model, controller and forecaster 
are all the same as the previous cases. Figure 4.11 shows 
the process response under IMC with no SSF. Figures 4,12 
and 4.13 show the cases where the linear forecaster and the 
AR forecaster are used respectively. A quick comparison of 
the three figures seems to show that there is no obvious 
advantage in using SSF for noisy signals and that the input 
Sequence has become more oscillatory in cases with SSF. 


However, a careful comparison of the output signal in 
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Figures 4.11 and 4.13 shows that the output in Figure 4.13 
is randomly distributed around zero whereas the output in 
Figure 4.11 has certain trends possibly caused by the 
disturbances. Therefore, the use of SSF helps the process 
to compensate for disturbance which has a trend even’ though 
it may be noisy. 

A common solution to the noisy measurement signals is 


to use a low pass filter. The equation is 


F(z) = ————_—_—_— (4.17) 


Simulation for the previous three cases are repeated 
with the additionscfea low pass filter (a = 078) @ein the 
System. The results are shown in Figures 4.14 to 4.16. The 
filter seems to do very little in improving the control 
performance but it does help to reduce the oscillation in 
the input sequence considerably. No attempt was made to 
optimize the value of the filter constant a. 

In summary, it has been demonstrated through simulation 
examples that SSF can be used to improve disturbance 
compensaton for IMC. The AR forecaster is better than _ the 
linear forecaster in an environment with noisy signals. A 
low pass filter should be used to filter the noisy signal 


before is is used by the forecaster. 
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4.5 Conclusions 


1) 


2) 


3) 


4) 


In the feedback interpretation of the IMC structure, the 
design rule that the IMC controller should approach the 
inverse of the process model is readily seen as being 
equivalent to the classical design observation that 
perfect control is approached as the feedback gain 
approaches infinity. The use of high gain for perfect 
CONEnOM Ineereedbackscontrol is srestrictedsbyestability 
and performance considerations. Similarly, the inverse 
model controller in IMC has the additional requirement 
of a perfect model in order to provide perfect control. 
IMC can be interpreted aS a classical feedforward 
control scheme based on the forecasted disturbance 
rather than the actual measured value. This 
interpretation helps to explain the function and form of 
the controller, the filter and the reference model. 

IMC can be interpreted as an SSF system that uses’ the 
current value of the residual as an estimate of the 
future value. The use of other forecasting methods’ can 
therefore be seen as an attempt to provide more accurate 
forecasts. 

Simulation examples show that the addition of linear or 
AR SSF to IMC improves its ability to compensate for 
unmeasured disturbances. However, the advantage gained 
is that the disturbances will be compensated A sampling 
periods earlier where A equals the sampling plus process 


delays. Moreover, the improvement is easily destroyed 
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by stochastic noise which introduces uncertainty into 
the forecasting. Therefore, the practical justification 


of using SSF with IMC is limited. 


5. SSF and the Adaptive Predictive Control System (APCS) 


5.1 Introduction 

The APCS is a multivariable adaptive control system for 
the control of linear, time-invariant processes with unknown 
parameters and subject to "bounded" stochastic noise and 
unmeasured disturbances. The basic theory of APCS was 
published by Martin-Sanchez in 1976 based on the work of his 
doctoral. thesis in 1974. An extension of the results to 
handle multivariable processes with time-delays was included 
in an US patent on APCS filed by Martin-Sanchez in 1977. In 
these early works, the stability of APCS was based on 
Popov's hyperstability criterion (1963). A new proof of 
convergence and global stability for APCS was reported by 
Martin-Sanchez, Shah and Fisher (1984). The new result 
involves the modification of the original adaptive mechanism 
to include a criterion for turning parameter identification 
on and off. This result has been further extended to cover 
time-delayed systems by Martin-Sanchez (1984). 

The strength of an adaptive controller comes from its 
ability to adapt its parameters continuously so as to 
improve control performance. The feature of stopping 
identification in APCS, though providing some distinct 
advantages, does raise doubts since the advantages of 
adaption are absent when the identification is off. This 
will be of particular concern if the criterion of stopping 


identification is mot chosen appropriately. The incentive 
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to improve the control performance of APCS was the original 
motivation for the work of this thesis. The solution is to 
make use of the estimation error or residual in the control 
Calculation. 

The purpose of this chapter is to investigate the use 
of SSF on the estimation error in APCS for compensating 
unmeasured disturbances and/or modelling error during the 
period when APCS parameter identification iS suOUt. 
Simulation examples will be used as illustrations for the 


ideas presented in this chapter. 


5.2 Description of APCS 

The Adaptive Predictive Control System is a globally 
Stable control system for multivariable, stable-inverse, 
time-invariant processes with pure time delays and in the 
presence of bounded unmeasured disturbances plus process and 
measurement noise. The following is a brief description of 
the control system for the case of Single-input 
Single-output processes. Fig. 5.1 shows a schematic diagram 


of the system to be described. 
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Consider the single-input single-output process in 
Figure 5.1. Assuming there is no measured disturbances, the 


process can be represented by a discrete ARMA model 


Vk) Papo oOy Gkedad)mts 640, Ck-d— imei (ik) C5205) 


where 


Sy (kd year (eave yee Ke O20 meter 
Cini a) WORF ong ) (552) 


@,(k-d-1) is a vector of past actual values of the 
Process iMoutputy Lyq,;imand)Pthek=control*® anputsh) fo. The 
dimension of ®, depends on the assumed process model order. 
ad is an integer which represents the pure time delay of the 
process. £&(k) represents the effect of unmeasured 
disturbances on the process output at the kth sampling 
instant. 9,5 and @; are unknown vectors of process parameters 
of appropriate dimension. 

Because of measurement errors and noise, the measured 
process variables are different from the actual values. 


Therefore, 


y (k) aa Se deg OS) Wo) 


u(k) = u,(k) + nylk) (onay 


(1.8) 


(8.2) 
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The corresponding measured ® vector becomes 


&(k) = @,(k) + ng(k) (5.5) 


Substitution of equations (5.3) to (5.5) into equation 


(5.1) gives the final form of the process model. 


y(k) = 6 x(k-d-1) + A(k) (5.6) 


where 


6 = [ 6) 6, ] = process parameter vector 
Reak—d—11)) == (SOO k-d-1)euCk-d=1)) J @= processm1, 0: vector 
Nee kK sO=l meee leng. (Kodo) niek>sO-i)) e=enovusemvector 


ACK) = on, (k)e— 62nehk—d-1) eet Ck) R=Bperturbacion term 


The one step ahead prediction of the process’ output 
¥(k|k-1), based on the estimation at time k-1 of the process 


parameters vector 6(k-1), is given by 


9(k|k-1) = 6(k-1) x(k-d-1) (5.7) 


The corresponding prediction error or estimation error 


is 


e(k|k-1) = y(k) - $(k|k-1) (5.8) 
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The estimated parameter vector is updated by the 
following recursive algorithm at every sampling instant 


depending on the value of w(k). 


Wilk) x(ked=1jetkikot) 
RS) SGM G Ry Ge ee Se (5.9) 
i+w? (k)x*t (k-d-1)x(k-d-1) 


The scalar wW(k) provides a basis for turning parameter 
adaption on and off. The definition of wW(k) is given as 
follows. 
avviak) ==" Omrieandsonly af 

HG) Kore lmesee Awa Ute A AK) Smec Agee (5.0109 


where 


2+2Wi2x*t (k-d-1)x(k-d-1) 
Ag'(W1,44,k) = ——————__ A, C5) 
2twi2xt (k-d-1)x(k-d-1) 


Wists lige Geta cence Oo eA eet A eee Xam) OK) y| 
0<k<@ 


A, is an estimate of a constant upper bound on the 
absolute value of the perturbation term A(k) for all k. 
Amn is the minimum value of this upper bound. 

by wit <) 2 (kK) Seek) askbme = @itandtonly if 
le(kiik-1)| > Ag’ (Wi,4a,k) 2044 (5a 25) 
where wWa(k) is defined as follows: 
i )euge(k) =ovpeeit. | etk}ksn)heseAgh wr, Agsk) 5h se 
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2+2W,,7x* (k-d-1)x(k-d-1) 
Se we AS —<—————————— (5.14) 
2+Wn2x* (k-d-1)x(k-d-1) 


i ig) 
2(|e(k]|k-1)|-Ag) 


Wa? (k) (S Als) 


TALE ERA ee ae) 


DESAI A ay ke) <r 6 (ki i <a re AG ke) 


The calculation of w(k) requires an operator-specified 
parameter Ag which is defined above. The definition of w(k) 
is closely related to the proof of stability which is given 
by Martin-Sanchez et al (1984). 

The prediction at time k of the process output at time 


Rec ery. (kts Kms sO 1Veraby 
7(k+dt+1/k) = O(k)x(k) = Oo(k)O(k) + 6, (k)u(k) fois 1) 


The control input u(k) is computed to make the 
predicted output y(k+td+1|k) equal to the desired output 


ya(kt+td+1). Therefore, 
(ki) = OU Klay a RO 1s) mee oak) oka ad (a5 NP) 


In order to guarantee a finite control input u(k), 
0,(k) must be non-zero. Martin-Sanchez et al (1984) have 


shown that 96,(k) can be made non-zero by an appropriate 


choice of w(k). 
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The control error is defined as 


e(kd@aty (kh -ay otk) (52480 


Equations @(5tijiteto (5918)8 desonibes@ag basic APCS 
algorithm. En@yEigees5Son ,m@etheggeneral formulation of ‘APCS 
also includes a "driver block". It can be interpreted as an 
extension to the traditional concept of the "reference 
model", At each sampling instant, the driver block 
generates from the operator-specified setpoint y, a desired 
HEOCeSGMOULDUSY GmmLOLea) LUtuLes samplingssinstant. This 
desired process output y,y belongs to a desired process 
output trajectory that satisfies a specified performance 
criterion. The driver block provides a means to enable the 
control algorithm to handle non-minimum phase systems and 
variable time-delayed process. 

As a summary, the characteristics of the APCS are 
listed as follows. 

1) It uses an ARMA process model which is identified 
on-line. 

2) It makes use of an I/O process model to calculate 
control action such that the predicted output is equal 
to the desired output. 

3) It has a criterion for stopping parameter adaption. 
This requires an operator-specified parameter which is 
an estimate of the upper bound of the perturbation term. 


4) Global stability is established for multivariable, 
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Stable-inverse, time-invariant processes with pure time 
delays and subject to bounded disturbances. 

5) It contains a driver block which generates a desired 
output from an operator-specified setpoint, the current 
parameter estimates, and in some cases an objective 


function and/or design method. 


5.3 Discussion of APCS 

The basic control algorithm of APCS is quite similar to 
the Self Tuning Regulator (STR) proposed by Astrom et al 
(1973)% However, it possesses several advantages over the 
STR. The first advantage of APCS is the guarantee of global 
Stability. Since the identification algorithm in most 
adaptive control systems is a highly nonlinear algorithm, 
the BproofetofimstabilityGis veryeddifeculty ERoregsSTRyaiu can 
only be shown that the algorithm will converge in most 
cases. In contrast, global stability has been proved for 
APCS. 

The second advantage of APCS is the criterion for 
turning parameter identification on and off. This not only 
saves computation time, but also avoids problem like 
estimators windupeewhichs fiSeecommon sinegescRs The ehird 


advantage of APCS is that the assumption it makes about the 
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unknown disturbances is very general. Most of the other 
adaptive control algorithms impose more restrictive 
assumptions on the disturbance or noise dynamics. 

The APCS also has some weaknesses. The following 
discussion will focus on two problems associated with the 
APCS adaptive algorithm and the feature of stopping 
parameter identification. Implications from these problems 


and solutions to them will also be discussed. 


5.3.1 APCS Adaptive Algorithm 

The APCS adaptive algorithm is a Dawa hy, simple 
algorithm to implement and the prediction ¥(k+A|k) 
calculated using the estimated model converges rapidly to 
the true value. However, the convergence of the model 
Parameters can be quite slow. This is due to the form of 
the adaptive algorithm (equation (5.9)) as explained by Song 
(1983). All the parameters adapt with the same gain which 
is the estimation error. The direction for adaption for 
each parameter depends on the sign of the elements in _ the 
vector x(k-d-1). This is different from the recursive least 
square algorithm which allows the parameters to adapt in 
different directions and with different gains because of the 
covariance matrix. This slow parameter convergence is not a 
problem if» the model is used for prediction since accurate 
prediction does not require accurate parameters. However, 


it becomes a problem when the process model is used to 
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calculate control action. In equation (5.11), the parameter 
6; acts like an overall gain and affects directly the value 
C@pathercontrolvachi ongcaliculated. sett 64 15 smalleand close 


to zero, even a slight change in its value will affect the 


Gon tsoli¢act vonmesigni fircant ly: This may lead to erratic 
controlkmaction twhich “is? undesirable in most cases. 
Thenvekorey, Sthegeaccuracym of PtGqi bei setverymnimpontant a One 


solution @tomresolve Thhiswepracticaheepnoblem@nisUetoetrirst 
identify 0, accurately off-line and then use it as a known 


parameter in the process model. 


5.3.2 On/Off Feature for Identification 

An adaptive control algorithm works by constantly 
adapting the parameters of the system so as to minimize the 
error between the output and the setpoint. When the 
parameters are not being adapted, the controller is fixed 
and it loses its ability to improve its performance at a 
certain operating condition or to maintain the level of 
performance when the operating condition changes. 

In APCS, the parameters of the process model are being 
adapted until the estimation error is smaller than A,'. A,' 
is a function of A,, which is an estimate of the upper bound 
of the absolute value of the perturbation term. In essence, 
APCS will stop parameter identification when it is not sure 
whether the eStimation error is due to modelling error or 


disturbances. This feature provides the advantage that the 
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model parameters will not be corrupted by the disturbances 
and the model will represent the true input-output 
relationship of the process. 

However, since the control error is equal to the 
estimation error, it means that APCS will only try to bring 
the output to within A,' of the setpoint. Therefore, the 
choice of A, dictates the ultimate performance level of the 
control algorithm. In practice, since the magnitude of the 
disturbances is unknown, A, has to be estimated. In order 
to guarantee stability, the choice of A, has to be lenient. 
This results in degradation of the control performance. 

The estimation error contains information about’ the 
unmeasured disturbances and modelling error. When the 
identification is on, this information is utilized since the 
estimation error is used in the adaptive algorithm to update 
the I/O model parameters. When identification is off, the 
estimation error is not used and the information it contains 
is wasted. Therefore, an obvious improvement to APCS is. to 
make use of the estimation error for control when the 
identification is off. Two cases will be described in the 
following paragraghs to show potential improvement through 
using the estimation error for control when the 
identification is off. 

Consider the first case where a process is subject to 
some unmeasured disturbances and a non-zero A, is chosen for 
the identification algorithm. Assume that when the 


identification stops, the model is perfect and the 
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estimation error is equal to the disturbances. Since the 
estimation error is not used when the identification is Ofte 
information about the disturbances is only fed back to the 
controJler through the use of measured process output 
values. Compensation of the disturbances is slow and 
requiresethatm the controller ycontainse inherent esintegral 
action. Alternately, by feeding the estimation error back 
to, the controller; “more direct’ control -action» ‘can be 
generated and the disturbances can be compensated faster. 
Consider the second case where a non-zero A, is used 
but for a while there is no disturbance affecting the 
system. Identification stops when the estimation error is 
Smaller than A,g' but the process model is not perfect. 
Theoretically, identification should be allowed to continue. 
However, Ag cannot be made zero because the disturbance is 
unknown. Therefore, the modelling error is left 
uncorrected. An offset will be present in the output. By 
making use of the estimation error or residual in the 
control algorithm to adjust the control action, the offset 
can be eliminated even though the modelling error is not 


corrected. 
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5.4 SSF and APCS 


5.4.1 Description 

Tt has been explained in the previous section that the 
control performance of APCS can be improved, when the 
TACTIC LciCcatiOnmeusee off, | by =makingususemson thesestimation 
error. This section will describe how the estimation error 
can be used and how SSF fits in this additional feature for 
APCS. 

Depending on the accuracy of the process model, the 
estimation error contains unmeasured disturbances and/or 
modelling error. This error is already present in the 
process output. If the value of the error is used directly 
by the controller, the control action, whose effect will be 
felt some time in the future, will be correcting the current 
error at a future time. What is really needed is the value 
of the future error so that the effect of the current 
control action can be chosen to minimize with the future 
error. Therefore, instead of feeding the estimation error 
directlytback itoathemscontrollempairts presenttand past values 
are, usedre:in SSSEwitomi forecast Gay future errors for sthe 
controller. 

Recalling the feedforward control system based on 
estimated disturbance in Chapter Three, the situation is 
Similar here for APCS. SSF can be used to forecast future 
values of the eStimation error based on the current and past 


values. The general form of the forecasting equation is 
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e(ktdt+1) =i f (e(k), e(k-1) 220) Cs TR) 


The forecasted value can be used by the APCS controller 
which is already an inverse model controller. Therefore, 


equation (5.17) is modified to 


u(k) = 6,°'(k) Lyg(k+d+1) - 6(k+d+1) - 05(k)@(k)](5.20) 


The modified APCS structure with SSF is shown in Figure 
5.2. The forecaster is designed to be used when the APCS 
identification is off. When the APCS identification is on, 
the forecaster has to be turned off so that it will not 
affect the identification. In order to provide a "bumpless" 
on/off switching of the forecaster, the forecasted value is 
multiplied by a switching factor to switch the forecaster on 
and off gradually. The switching factor w which lies in the 
range of 0 to 1, is defined by the following equations. 
SWitchingpone mor 08+ey( 1a) C52 2.18) 
Switching off : w = (1-7)w Gar 2iib) 
where y is a parameter which determines’ the 
LatemOLeGwitGhingeand 0m sey ass 
y = 0 means no switching action 
y = 1 means bang-bang switching 
With the switching feature for SSF, the identification 
and control action of the basic APCS is not affected. 
Therefore, the stability proof of APCS is maintained and the 


modified APCS (with SSF) is also stable. 
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Before proceeding with some simulation examples of SSF, 
an alternative to the second level of control by SSF will be 
discussed. The alternative is to allow parameter 
Dene icicat’ OnueLOmerinwateal atime =—aTniseapproacheismfound 
in most adaptive control schemes which do not have on/off 
Switching for parameter identification. 

The disadvantage of this alternative is that the 
process model parameters could be corrupted if the residual 
error is due mainly to disturbances. Consider the situation 
where an accurate set of model parameters have been 
obtained. If identification is allowed to continue and some 
new disturbances come in, the parameters will be corrupted 
while the disturbances are being compensated for. Pia 
Setpoint change is required at this point, the response will 
be poorer with the corrupted set of parameters. Moreover, 
if this alternative is adopted for APCS, the conditions 
required fore@ ther stability sprooty will, “notes bem mete and 
Stability is not guaranteed. 

Therefore, the use of SSF to provide a second level of 
control can reduce the control error and prevent unnecessary 
COLLUDL ION OL the = DLOCceSssS 1/0 model by unmeasured 


disturbances. 


5.4.2 Illustrative Examples 
Some simulation examples will be given to illustrate 


the use of SSF in APCS and to explore some of the properties 


fic\\ne «sued don of ds 
9 
sit ted: ai Suiiaci 
Leuwhizss afd 2% Betqi 
nottewste sit sahinnop i il 2 
need sired 21S ISme Teg, {shew do Bs 
an68 Sia ours sme bosotls ee anys 4 
Setquites <d Edbw wpe tams 6 aas (ak pac 22 . 
a 2: 103 bedgenp@go> sited... esaapeeneb ere 0 
[liv sqnoges3 sat ytndeg etd? te bot kype? ai. | 
yuevesscM  -evetiekag fo 19@ besquaties ads @ 
anviziincy ois (AR vo ‘bengees | 2t- fortsencea 
bas son 06. so Cee doen —— of 


mae 
te, es on ee Wt 


ies year oe eA 307330 toxtnos ed: a2uhs?: ; 
oe a2 o 
betueseme «= yo febem O\l sepoeag. ons 


sxetzeulit 09 nevi 90 tthe selamexs nois ake 
srirseauia oid 3 enor sr0tgxe 62 bee 2084 nt 22 
7 : _ 7 7 : we 


Oy 


of the system. The second order process which has been used 


in 


Chapter Four will be used again. The process can be 


described by the following ARMA equation. 


y(k) = 1.030) y(k=1)e= 08 2898y(k=2) 
+ 0.158 u(k=1) + 0.101 u(k-2) (5s 22) 


A second order ARMA model as described by the following 


equation is used by APCS. 


y(k) = diy (koi tae Vik 2) tbe kode beck 2) (Se 2380 


The following cases have been examined via simulation 


studies. 


1) 


2) 


In the first case, the AP model has perfect parameters 
and no identification takes place. The process is 
Subject to a noise-free disturbance sequence. This is 
the first case described in Section 5.3.2. Figure 5.3 
shows the response of the system under APCS control. In 
Figures 5.4 and 5.5, a linear forecaster and an AR 
forecaster were added respectively. The figures show 
that APCS performs poorly in compensating for unmeasured 
disturbances and the addition of SSF improves’ the 
control performance significantly. 

In the second case, the AP model parameters are 
considered to have converged under certain operating 
conditions but they are not equal to the true 


parameters. A setpoint change is introduced. The A, is 


tfs 2) 


soijeivere 


ei 27399075 
ai etat 
2 enopis . : 
| > ea ee 
ht ,ttithod 6On: ‘Sibeat on ste | Spal: 
“4 ns 6ns sadeaietes taondit J 
wore zysupht eit smog 
dbercaseimy 102 pakasanoqnos ar 


3) 


98 


non-zero but there is no unmeasured disturbance. Figure 
5.6 shows the response of the system under APCS control. 
A steady state offset is present because of the 
imperfect parameters and a non-zero A,. Rigucemmo. 
shows the addition of SSF helps to remove the steady 
state offset. 

In real applications, the model parameters are generally 
not equal to the true parameters and there are 
unmeasured disturbances. In this case, the simulation 
was done with an imperfect AP model to compare the 
performance of APCS with and without SSF for regulatory 
Control. A noise-free disturbance sequence is used for 
Figures 5.8 and 5.9 while a noisy disturbance sequence 
is used for Figures 5.10 and 5.11. The results show 
that the modelling error does degrade the control 
performance of “APCS: The addition of AR SSF helps to 
compensate for both disturbance and modelling error. 
Moreover, the noise has no significant effect on the 


performance of the AR SSF. 
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5.5 Conclusions 


1) 


3) 


4) 


Adaptive control systems, such as _ APCS, use the 
Magnitude of the estimation error as a criterion for 
Switching identification on and off to prevent the 
disturbances from corrupting the identificaton of the 
true input-output relationship of the process. However, 
this criterion also has the disadvantage of not 
utilizing the information in the estimation error when 
the identification is off. 

The addition of SSF to adaptive control systems such as 
APCS provides a second level of control when 
identification is suspended. SSF generates future 
values of the residual for the controller to provide 
additional control action to compensate for unmeasured 
disturbances and/or modelling errors. The second level 
of control also helps to prevent unnecessary corruption 
of the process I/O model. 

When SSF is used with APCS, it provides the advantage 
that the stability proof of APCS is maintained since the 
forecasting is done when APCS identification is off. 
Simulation results confirm the improvement expected from 
the second level of control based on SSF. Modelling 
errors degrade the performance of APCS but can be 
compensated for by SSF. The performance of the AR 
forecaster is not significantly affected by the presence 


of noise and is therefore recommended for use over other 


forecasters. 
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6. Application of SSF 


6.1 Introduction 

In the previous chapters, the idea of using SSF as a 
means to compensate for unmeasured disturbances has been 
introduced. Numerical examples were used to illustrate the 
concepts presented. In thisfehapter, SSE is applied to a 
Simulated double effect evaporator in order to study its 
practical usefullness in a more realistic application. SSF 
is applied as an improvement to APCS for the compensation of 
unmeaSured disturbances. The performance is compared to 
that of the conventional feedback control and the basic 


APCS. 


6.2 The Evaporator 


6.2.1 Description of the Equipment 

The evaporator whose model is used in the simulation 
studies is the double effect pilot plant evaporator in the 
Department of Chemical Engineering, University of Alberta. 
The evporator has been described in detail by Fisher and 
Seborg (1976) and its schematic diagram is shown in 
Figure: 6.51% 

The first effect of the evaporator is a natural 


circulation calandria type unit with thirty-two 18 inch 
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long, 3/4 inch OD tubes. It is heated by process steam. 
The second effect is a forced circulation long tube unit 
with three 6 feet long, 1 inch OD tubes. It is operated 
under vacuum and is heated by the overhead vapor from the 
first effect. 

The evaporator is fully equipped with industrial 
electronic instrumentation and is interfaced to an HP/1000 
computer via a LSI-11 microcomputer in the Data Acquisition, 
Control and Simulation (DACS) Centre of the Department of 
Chemical Engineering at the University of Alberta. 
Therefore, evaluation of advanced control techniques can be 
eaSily carried out on the evaporator. 

The primary controlled variable of the evaporator is 
the product concentration, C2, from the second effect. lee 
is controlled by manipulating the steam flow rate, S. Other 
important controlled variables are the first effect holdup, 
W,, and the second effect holdup, Wz. They are controlled 
byemanipulating thes firsts effect e@boctom flow race oy seeand 
the second effect bottom flow rate, Bz, respectively. When 
the evaporator is used as a_e single-input-single-output 
system, the C2/S loop is used and the other loops are closed 


by conventional feedback control. 


6.2.2 Evaporator Model 
Several models have been developed for the double 


effect evporator in previous studies Newell (1971); Wilson 
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(1974). They range from a tenth-order nonlinear model to a 
first order transfer function model. The evporator model 
USed any this work @as the) fifth corder nonlinear model 
developed by Newell (1971). The model consists of five 
Ordinary differential equations and a set of algebraic 


equations. The equations are given as follows. 
Steam chest ;: 
SH Un AG Tis) (6.1) 


Brvrst) effect 


aw, 
So Se PY Be ratOn (6.2) 
dt 
dc, 
W; ——- =F (Gn - Cw rT 0,C, (6.3) 
at 
dh, 
W, —— =F (ho-h,) - 0,(Hv,-h,) + Q, - L, (6:..4)) 
ohn 
where 
Oy =2.U ACh) oR) 
Ome (Oo +o) ea Vance) (6.6) 
Oona Up Avon Toe) (36 wath) 
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Second effect : 


aw 2 
Gel eshte Ueley St ON (6.8) 
dt 
dGe 
We = BaiG. Gs) + 02C>2 C659) 
at 
where 
dhe 
O27 (Hv; -ne= Ga) =F O02 ethane 
0C>2 
dhe 
ar By(Cs-G,) (6 10) 
0C2 
Some property relations are 
Hv = 0.4Tv + 1066.0 (6s) 
ite i Gh = 2.) UOwa2) 
HSOle=aTSOlelm—" (0. (6CSOL) ma ma2e COs!) 


A FORTRAN evaporator simulation program was written 
based on the above equations. The differential equations 
were solved by the Fourth order Runge-Kutta Method with a 
Simulation interval of 0.1 minute. Steady state operating 
data for the evaporator were obtained from Wilson's thesis 


(1974) and are tabulated in Appendix A. 
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6.3 Practical Considerations 

In applying APCS with SSF to a real process or a model 
of one, there are many practical considerations in the 
implementation. The following discussion will address these 
considerations in the application of APCS with SSF to the 


double effect evaporator model. 


1) Process model order 

APCS requires an ARMA process model whose parameters 
are identified on-line. The order of the model _ should 
theoretically match the order of the process. However, in 
practice, a low order model is uSually sufficient as an 
approximation. For the evaporator, a second order model is 


used. The equation of the model is 
Vik) me an VekKetie tide y UK—>2) ete by uk 10" fabou GK~2 me Gea) 


There is considerable modelling error as can be seen by 
comparing equation (6.14) versus equations (6.1) through 
(6.13). Note that one of the purposes for using SSF is to 


help compensate for modelling errors. 


2) Initial parameters 

The choice of the initial parameters for the model is 
critical to the performance of the control algorithm. A 
discussion on the choice of initial parameters for APCS has 


been given by Song (1983). In the following application, a 
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set of identified parameters from Song's doctoral thesis 
(1983) were chosen. They were used as starting values ina 
period of parameter identification which yielded a more 
accurate set of parameters for use as initial parameters in 


the simulation studies. 


3) Scaling of variables 

The evaporator model generates values in absolute 
engineering units. If the variables used in the process 
model are in engineering units, the magnitude of the model 
parameters may vary greatly and computational problems may 
result. Therefore, the input and output variables have to 
be normalized. In this case, the normalization is based on 


the steady state values. 


4) Sampling interval 

The choice of sampling interval for adaptive control 
systems should relate to the time constant of the process. 
For the evaporator, the sampling interval should be about 3 
to 4 minutes. However, the abies evaporator is usually run 
withealosampling! interval® ofem64 seconds: Therefore, a 
samplingawintervale@mot (ei) *minute is,used for the simulation 


studies. 


5) AR forecaster order 
Theoretically, the order of the AR forecaster should 


match the order of the residual series. However, since only 
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short term forecasts are made, a third order AR forecaster 
is used to reduce the computational effort. The form of the 


AR forecaster is 


BCk bi)ieata pol Ket kar Chait) page Gk 2) CGu15) 


6) Initial parameters for the AR forecaster 

The choice of the initial parameters for the AR 
forecaster affects the accuracy of the forecast. Using 
engineering judgment, the initial parameters have been 
chosen so that the AR forecaster is initialized to use the 
current value as its forecast. This is a good approximation 


and provides reliability for the SSF. The parameters are 


7) Identificaton algorithm for the AR forecaster 

Many parameter identification algorithms are available 
for on-line identification. The APCS projection algorithm 
(equation (2.9)) has been chosen because of its simplicity 


and the author's familiarity with it. 
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6.4 Simulation Results 

The performance of conventional feedback control, APCS 
and APCS with SSF was studied. The evaporator was subject 
to a noisy disturbance sequence in the feed flow rate. The 
control objective for the system was to maintain the product 
concentration at the setpoint value. 

Figure 6.2 shows the case for a well tuned feedback PI 
GOntrobler g\Kcou=—7200, K;, = 001). ‘The control action us not 
sufficient to keep the output at the setpoint. The product 
concentration deviates from the setpoint due to the feed 
disturbance. Figure 6.3 shows the performance of a well 
EUneC mm TeecabackmariD =controller (Kc = —S00sweK i = 105 157 
Ky = 0.3). The addition of the derivative mode improves 
control Since it provides anticipatory action (c.f. Chapter 
Three). However, the product concentration still deviates 
from the setpoint slightly. 

Figure 6.4 shows the performance of APCS. Since the 
identification is off, the controller parameters are fixed 
and the controller cannot compensate for the varying 
disturbances. Figure 6.5 shows the case where SSF has been 
added to APCS. SSF provides additional control action when 
APCS identification is OLE The disturbances’ were 
compensated and there is significant improvement in the 
OuULDUL, TeSponse: 

Comparing the four cases in Figures 6.2 to 6.5, it can 
be concluded that the concept of SSF, whether applied in 


conjunction with a predictive: scontroller (APCS) or 97a 
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non-predictive controller (feedback control), can be used to 
improve a control system's performance in disturbance 
compensation and/or in accommodating modelling errors. 
However, the practical justification will vary in individual 
cases and the advantages gained will have to be evaluated 


relative to cost and overall system performance. 
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7. Conclusions 


7.1 Conclusions 


Thesrconclusioens® resultingetfromm™ ehiss ’studyeacaneibe 


divided into two categories: 


Al) 


A2) 


A3) 


A) those relating to the performance of SSF, and 
B) insights gained by comparing other classical and 


modern control techniques with SSF. 


SSFecaniresult intperfecticontrollificthegrprocess 617/70 
model and the forecast are perfect. However, the 
quality of control deteriorates due to factors such as 
large amounts of random noise, strong nonlinearities, 
which make forecasting difficult. 

The forecaster should be able to predict the future 
values of the residual time series A sampling periods 
into the future where A is the sum of the sampling plus 
process delays. Autoregressive forecasters of moderate 
order provided better performance in the simulated 
studies than Simple linear extrapolation or the 
asSumptionethatem kt Ajeet Uk )e 

Since a process I/0 model is_ required for the 
implementation of SSF and since a model based (e.g. 
predictive) ipcontrollerSygives;#)ink egeneral, better 
performance than, say, a PID controller, the base case 


for comparing the advantages of SSF should assume _ the 
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A4) 


A5) 


A6) 


B1) 


USemoOf Fagepredictivercontrol ler. 

Several simulated runs showed that SSF helped to 
compensate for unmeasured disturbances and/or modelling 
errors, However, when used with a perfect predictive 
controller, gthe Bbest lathatgiSShieucan “do@sristiprovide 
compensation for the effect of a disturbance at time k 
rather than A _ periods later. Thus@e@itse practical 
usefulness is limited. 

The addition of SSF to adaptive control systems such as 
APCS, which uses the magnitude of the estimation error 
aS a criterion for switching identification on and off, 
provides a second level of controls: When 
identification is stopped to prevent the I/O model from 
being corrupted by the disturbances, the second level 
of control comes in to provide short term compensation 
for unmeasured disturbances and/or modelling errors 
based on information of the residual. 

When SSF is used with APCS, it provides the advantage 
that the stability proof of APCS is maintained since 
the forecasting is done when APCS identification is 


Offs 


Classical proportional plus derivative feedback control 
can be interpreted as SSF of y(k) in conjunction with a 


proportional controller. 
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From the equivalence relationship of the feedback 
control structure and the IMC structure, the design 
rule that the IMC controller should approach the 
inverse of the process model is readily seen as _ being 
equivalent to the classical design observation that 
perfect control iS approached as the feedback gain 
approaches infinity. (There are of course stability 
and performance’ limitations with most practical 
systems.) 

IMC can be interpreted as an SSF system that uses the 
current value of the residual r(k) as an estimate of 
the future value r(k+aA). This implies that IMC can 
also be interpreted as a classical feedforward scheme 
based on the forecasted (estimated) disturbance rather 


than the actual meaSured value. 


7.2 Recommendations 


i) 


The primary concern in this thesis is to investigate 
the idea of using SSF to improve control. The issues 
involved in the implementation of SSF, such as_ the 
choice of process model, forecaster and controller, 
have not been addressed in detail and are proposed for 


future work. 


q¢sri idssé 21965 x6 Farle 
feaijnbiq = Jean. Aste 


4 Pie. 
ai? esau Jeet Merage peer Se 
to Stemives: GetVas Gade <abitee Wi Xo sulev dnetws 
ied MP sade a ak 2a -. tiie euhey eresul -ed3 — 


an se se brew ideeeet | jovi gaala one Besgnqzemi sc oats 
vadiet Sonsd SEER Nae aR 06s no beeed 


2s 


2) 


3) 


4) 
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Extension of the use of SSF in multi-input multi-output 
(MIMO) systems is recommended for investigation. 
Experimental evaluation of SSF should provide results 
BOEMEVerificatione¥of or the conclusions drawn from 
Simulations in this thesis. 

The interpretation of the forecaster as a disturbance 
model suggests a possible relationship with the 
Internal Model Principle by Wonham. Future work can 


also be done in this area. 
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Appendix A : Steady State Operating Data for the Evaporator 


The following steady state operating data for the 


double effect evaporator were taken from Wilson's thesis 


(1974)°and listed in SI units. 
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Twi 


Wi 


ei 


H1 


Tw2 


Bl 


B2 


CF 


HF 


steam temperature 

First effect tube wall temperature 
First effect holdup 

First effect concentration 


First effect enthalpy 


Second effect tube wall temperature 
Second effect holdup 

Second effect concentration 

Second effect enthalpy 


Condenser tube wall temperature 


Steam flow 
First effect bottom flow 


Second effect bottom flow 


Feed flow 


Feed concentration 


Feed enthalpy 
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20.64 kg 
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OS015S i kg/7s 
0.0263 kg/s 
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Appendix B : Data for simulation examples of Chapter 3 


General data 

Process : “y(k)) =)0s6y (kod) h+40. Gutkadjucee Ck) @+O0 S5eGka1) 
Constant setpoint sf y ,@=95.0 

Simulation time : 100 sampling periods 

Process subject to random disturbances 

Predicbive controllerg: sulk) m= sly. Akad )or (kt 1)-0.8y (kK))170.6 
binear SSF :2 (kl) -serek igre cUk ) eo Bo( kd) ) 


ARMSSRe sone (kta) B= tay ckhi)iet tagr Geei jt a5 (k-2) 


Initralevaluest)) a; (=u 
a2 = 030 
a3 = 0.0 


Identification : APCS adaptive algorithm 


Cases and specific data : 
1. Feedback controller 
Figure 3.6 Proportional controller 
Ke S017 
Figure 3. / Proportional-derivative controller 


Ke a ety, Ky cs ‘lee 0 


2. Predictive controller with and without SSF 
Bigune momo Perfect SSF 
FaGunees ac No SSF 
Figure 3.10 Linear SSF 


Figure: 3.1 |) eARPos® 
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Appendix C : Data for simulation examples of Chapter 4 


General data 
Process : y(k) = 0.158u(k-1) + 0.264u(k-2) + 0.226u(k-3) 
ae Oe lO6U (kK =4) et 0096UCK—5 eee 053u(k=6) 
FeO O27 Uk ar) eer Ou CK=o)) 
Constant setpoint : y, = 0.0 
Tuning parameters : P=N = 8 
Y= Yoo Sees = Yo = | 
B, = B2 = ... = Be = 0 
Simulation time : 100 sampling periods 


ta Nearaecosk) sae niktl)) = r(k)o + o(r (ke — r (k=4))) 


ARaSSHe:) ptkti)ie= ark) + asrik=t) = asr(k-2) 


Thittl sleeveless sn, Van. C 
icy i= ORO 
(Qf) 0.0 


Identification : APCS adaptive algorithm 


Cases and specific data : 

1. Compensation of noise-free disturbances 
Figure 4.7 IMC 
Figure 4.8 IMC with Perfect SSF 
Figure 4.9 IMC with Linear SSF 


Figure 4.10 IMC with AR SSF 
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Compensation of noisy disturbances 
Figure 4.11 IMC 
Figure 4.12 IMC with Linear SSF 


Figure 4.13 IMC with AR SSF 


Compensation of noisy disturbances with filter 


i= 7 0: 
Filter : F(z) = 


Figure 4.14 IMC 
Figure 4.15 IMC with Linear SSF 


Figure 4.16 IMC with AR SSF 


Appendix D : Data for simulation examples of Chapter 5 


General data : 


Process : y(k) = TOS CR y Gk—1) Ma O S26 9m. (Ko) 
Pere Oe Uh GK—1t) ers meiO tk cae) 
Constantesetpormnte: yi =9 120 


Simulation time : 


Linear SSF : £(k+1) 


AR SSF r(k+1) 


rACk 


a,c (k) 


Initial values 


Identification 


Cases and specific data 


Te 


model 

Figure 5.3 APCS 
Figure 5.4 APCS 
Riugures 5.0 APCS 


Setpoint tracking 
Setpoint change 
Figure 5.6 APCS 


Figure 5.7 APCS 
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of noise-free disturbances with imperfect 


APCS 


APCS with AR SSF 


of noisy disturbances with imperfect model 
APCS 


APCS with AR SSF 
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Appendix E : Data for simulation examples of Chapter 6 


General data 

Process : 5th order nonlinear evporator model 
Simulation time : 120 minutes 

Process subject to feed disturbances 


ARS S Kae: r(k+1) = a,r(k) ct aor(k-1) a QaOGk= 2) 


Tikstie evad UGS es) mos = 10 
($65) @.0 
(@ocy OF 0 


Identification : APCS adaptive algorithm 


Cases and specific data : 
1. Feedback control 


Figure 6.2 Proportional-integral controller 


Ke=200 Kk, =0..1 
Figure 6.3 Proportional-integral-derivative controller 
Ke=300 K,;=0.15 Ky=0.3 


2. APCS with and without SSF 
Figure 6.4 APCS 


Begu nemo. APCS with AR SSF 
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